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Abstract 

After a brief review of the derivation of the standard phase time formula, based on the use of 
the stationary phase method, we propose, in the opaque limit, an alternative method to calculate 
the phase time. The new formula for the phase time is in excellent agreement with the numerical 
simulations and shows that for wave packets whose upper limit of the momentum distribution is 
very close to the barrier height, the transit time is proportional to the barrier width. 



I. INTRODUCTION 

a ' 

The time spent by a particle to tunnel across a barrier surely represents one of the most intriguing 
and challenging discussions found in literature. After the stimulating articles of MacColl [T] and 
Hartman [2] on the dynamics of the wave packet which tunnel potential barriers, many tunneling time 
OA , definitions have been introduced and paradoxical effects, such as supcrluminal velocities, discussed. 

Of special interest for us is the discussion on the time spent by non-relativistic particles to cross the 
classical forbidden region. The extensive literature on tunneling times is reviewed in many reports. 
qq ■ For a detailed discussion on phase, dwell and Larmor times we refer the reader to the report of Hauge 

and Stovneng [3J, for photon and particle tunneling reviewed by a unified time analysis to the report 
of Olkhovsky, Recami and Jakiel (4], and, finally, for a clear, comprehensive and complete discussion 
of the tunneling time definitions, paradoxes and proposed solutions to the excellent work of Winful [5]. 
In this report, it is also found a challenging clctromagnetic analogy with the frustrated total internal 
reflection and resonant tunneling. 

In this paper, we present a detailed analytic and numerical analysis of the phase time for non 
rclativistic wave packets which tunnel opaque barriers. In the opaque limit, due to the filter effect, 
approximations on the transmission coefficient allow to find a closed formula for the time in which 
\ the peak of the transmitted wave appears in the free region after the barrier. The new formula, 

which generalizes the well known formula obtained by the stationary phase method, shows that, for 
momentum distributions whose upper limit is very close to the barrier height, the phase time is 
proportional to the barrier width. The study is done for a potential barrier V with discontinuity in 
x = and x = L. 

The method commonly used in calculating the transmitted amplitude is based on solving, sep- 
arately, the Schrddinger equation for stationary states within the potential region and in the free 
regions before and after the barrier and, then, imposing the continuity of the wave function and its 
derivative at the discontinuities of the potential. After simple algebraic computations, we find the 
following transmitted amplitude [5j[6] 



T(k) = 
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where q — v 'w 2 — k 2 and w = \J2mV a jh. The resultant transmitted wave packet is obtained by 
integrating over all the possible stationary states modulated by a weighting function g(k), 

rkM 

§ T (x,t) =N dkg(k) \T(k)\ e tip e nH*-L)-Et/h] ^ ^ 

J 

where tp~kL, with ip = arctan[ (2 k ~w 2 ) tanh(gL) / 2 k q], is the phase of the transmitted amplitude, 
k M = y/2 mE M /h is the upper limit of the momentum distribution, and TV is a normalization constant 
containing the information of the number of incoming electrons. The condition k M < w guarantees 
that the allowed energies are truncated before or at the barrier height V . We do not have above barrier 
contributions and, consequently, only tunneling is responsible for the transmitted wave. We restrict 
the discussion to pure tunneling, i.e. E M < V , to avoid the phenomenon of multiple diffusion [TIE] ■ 

As is well known [H[S], the use of the stationary phase method allows to calculate the time in 
which the peak of the transmitted wave appears in the free region after the barrier without explicitly 
solving the integral of Eq.([2]). Unfortunately, this is not sufficient for determining the transit time. 
In fact, the use of the standard phase time formula requires a careful analysis on the applicability of 
the stationary phase method. Without such an analysis its indiscriminate use could result in wrong 
theoretical interpretations on the dynamics of the particles tunneling. 

In the next section, we briefly revise the standard derivation of the phase time based on the use of 
the stationary phase method. Then, in section III, we obtain, by an analytic study of the transmitted 
wave, a new formula for the phase time . This is done in the opaque limit. The discrepancy between 
the standard and the new formula for the phase time is clear for E M close to V . To confirm the validity 
of the new formula, a comparison between analytical results and numerical data is presented in section 
IV. The agreement is excellent and suggests the use of the new approach proposed in this paper as 
a new method for estimating the phase time. Our conclusions and possible future investigations are 
drawn in the final section. 



II. REVISING THE STANDARD PHASE TIME FORMULA 



The old question of tunneling times is often addressed by studying the phase time through the use 
of the stationary phase method [HUH]- This method provides an approximate way to calculate the 
maximum of an integral. The main idea is that sinusoids with rapidly changing phases will add 
destructively. This basic principle of asymptotic analysis allows to find the maximum of an integral 
independent of the details integrand shape. Such a maximum depends on the derivative of the 
integrand phase calculated at the mean value of the wave number. For the transmitted wave given in 
Eq. ([2]) , the phase term is stationary for 



EC 



ip + k(x — L) — 

n 



= , (3) 

(k) T 



where the prime stands for the derivative with respect to k and 



rk M I rk M 

(k) T = dkkg 2 (k)\T(k)\ I dkg 2 (k)\T(k)\ . (4) 



By using Eq. (|3]) , we find that the phase of the transmitted wave is stationary at x = L for times 
which satisfy 

Et k [w 4 sinh(2 qL) + 2 k 2 (w 2 - 2 k 2 ) qL] ' 



H q [w 4 cosh(2 qL)) + 8 k 2 q 2 - w 4 } 



= . (5) 

(k) T 



For thin barriers (k M L 1) the modulus of the transmitted coefficient is close to 1. The trans- 
mitted wave packet has the same form of the incident packet and, consequently, (fc) T ~ /c , where ko 
is the center of the incoming momentum distribution g(k). The phase time is then proportional to L. 
Observe that for very thin barriers we can always guarantee q L <C 1. 
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For much thicker barriers (k M L 3> 1), we enter in the so-called opaque limit. In this limit, the peak 
of the transmitted momentum distribution is shifted to higher wave numbers. This effect is known 
in literature as filter effect Before beginning our discussion on the phase time formula ([5]), let us 
briefly discuss the filter effect. In the opaque limit, the modulus of the transmitted coefficient can be 
approximated by 

\T(k)\ « Akqe~ qL jw . (6) 

By using this approximation, which also implies g(k)\T\ ps g(k M )\T\, and changing in Eq. ((4]) the 
variable of integration from k to g, we obtain 



(k) T 



I dqk 2 q a e~ 2qL / / dqkq a e - 2qL 

J Q M J Q M 



(7) 



where q M = y/i 



' qM J qM 

k 2 M . The integrands can be now expanded in series around q M , 



k 2 q 3 



Kill, + (3» 2 - 5gi f )gL (q 







(q 



k q 3 = k M q 3 M + (3w 2 -4 q 2 M ) q 2 M (q - q M ) / k M + O (q - q M ) 

Observing that the main contribution to the integrals in Eq.([7|) comes from the lower limit q M , we 
obtain 

k 2 M q M + (3w 2 -5q 2 M )/2L q h 



(k) T 



k M q M + {3w 2 -4q 2 M )/2k M L 2k M L 
The phase of the transmitted wave is then stationary at x — L for times which satisfy 

Em tspM ^ k M 



(8) 



(9) 



In the next section, we shall propose a new method to calculate the phase time. The new phase time 
formula which reproduces Eq.((9]) for q M ~ k M foresees transit times which arc proportional to the 
barrier width L for q M k M . 



111. PROPOSING A NEW PHASE TIME FORMULA 

In the previous section, we have estimated the position of the maximum of the transmitted wave 
packet by using the stationary phase method. In this section, we propose an alternative method to 
calculate the phase time formula. For opaque barrier tunneling it is possible to calculate explicitly 
the derivative with respect to time of the electronic density at the barrier edge and finding when it 
is equal to zero. This allows us to obtain a new formula for the time in which the maximum of the 
transmitted wave packet is found at x = L. 

For opaque barrier tunneling, in solving the integral which appears in Eq.([5]) we can use the 
approximation given in Eq. (|6]) and change the variable of integration from k to q. Consequently, the 
expression of the transmitted wave at the edge of the barrier (x — L) becomes 

$ T (L, t) w 4 N g{k M ) j dq q 2 e~ qL e 4 (v>-Et/n) j w 2 _ ^ 

J qM 

Due to the filter effect, the phase ip and the energy E can be expanded as follows 
<P = <Pm — 2(g— q M ) I k M - q M (q- q M ) /k M + {q-q M ) 
E = E M - K 2 q M (q- q M ) / m- h 2 (q- q M f / 2m . 

Introducing the new adimensional variable p(q) = (q — q M )/k M and using the previous expansions, we 
obtain for the electronic density at x = L, 



|$ T 0M)l * 16|iV| k M g 2 (k M ) 



p(w) 



dp p 



exp{-p k M L + i [ a(t) p + fi(t)p 2 ] } 



S(t) 
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with 



a(t) 



Qm^m t 
h k M 



and P(t) 



E M t q h 



(11) 



The subject matter of this section will be the accurate analysis of S(t) and the calculation of its 
derivative with respect to time. A first approximation is to consider the first terms in the expansion 
of the phase time exponential, 



S(t) « | 8(0) + i [a(t) s(l) + P(t) 8 (2) ] - \ [a 2 (t) s(2) + 2 a(t)p(t) s(3) + f3 2 (t)s(4) 
where 



(12) 



s(n) 



dp p 



1* 



p n exp[-pk M L] 



(n + 2)\ 
(k M L) n - 



1 



(q M L) 



n + 2 (n + 2)(n+l) 



The case q A 



If we limit ourselves to the analysis of processes in which q M is of the order of k M , we find that j3(t) 
is of the same order of a(t). Observing that in the opaque limit s(n) 3> s(n + 1) and that the /3-term 
in S(t) is coupled to s(n) with higher n, we can approximate S(t) as follows 

S(t; q M - k M ) » | s(0) + % a(t) s(l) - \ a 2 {t) s(2) \ « s 2 (0) + a 2 (t) [ s 2 (l) - s(0)s(2) ] . 

Deriving g M ~ k M ) with respect to time and equating to zero, we obtain a(t) = 0. Thus, in this 
limit, we reproduce the well known stationary phase condition ((9|), 



St(t; q M ~ k M ) = 



= 



E M t 



(13) 



The case q M <C fc A 



In this limit, a ps —2 and ps V t/H. The /3-term cannot be neglected because for time of the order 
of KL/E M it becomes comparable to the a-term. This implies in our approximation that the terms 
s(n), ts(n + 1), and t 2 s(n + 2) are of the same order. Consequently, 



S(t;q M < fc M ) 



S (0) 



2 2 

4s( 2)-4 S (3)^ + S (4)^ 



_ 2s (l) + ^ s (2) 



2 2 

s 2 (0) + 4 [ S 2 (l) - S (0) S (2) ] + 4 [ S (0>(3) - S (1) S (2) ] ^ + [ S 2 (2) - S (0) S (4) ] 



By taking the derivative with respect to time and setting it equal to zero, we obtain 

V t 2 [s(l)s(2) - s(0)s(3)] 2wL 



S t (t;q M < k u ) = 



•(2)- a (0)*(4)] 
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(14) 



The transit velocity, defined as the ratio between the barrier width and the time in which the peak 
appears in the free region after the barrier, is then given by 



^tra 



9 Vo_ 
2 V 2m 



(15) 



This analytic result is confirmed by numerical calculations, see Table 1. The details of our numerical 
simulations are found in section IV. 
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The general case 



Observing that for increasing times the terms a(t)s(n), j3(t) s(n + 1), and /3 (t) s(n + 2) become 
comparable, we obtain for S(t) the following expression 

S(t) w s 2 (0)+a 2 {t) [ S 2 (l)-s(0)s(2)]+2a(t)/3(t) [ s(l)«(2)-«(0)*(3) ]+/? 2 (f) [ s 2 (2)-s(0).s(4) ] . (16) 

In deriving S(t), we use ctt{t) = 2q M E M /hk M and f3t(t) = E M /h, and after simple algebraic manipu- 
lations, we find 

g^w = 2 ( W /fc Af ) 2 [a(l)s(2) - a(0)s(3)] + 4 (<WM[s 2 (l) = «(0)a(2)] , s 

h [ S 2 (2)- S (0)s(4)]+4(g M /fc M )[ S (l) S (2)- S (0) S (3)]+4(g M /A :M ) 2 [s 2 (l)- S (0) S (2)] ■ 1 j 

For q M ~ fc M , remembering that in the opaque limit s(n) 3> s(n +1), we find 

gA, _ 4( gM /fc M )[5(l)s(2)-s(0)s(3)] = 
ft ~ 4(g M /fc M ) 2 [s 2 (l)- S (0)s(2)] fc M ' 

as anticipated by Eq. (|13[) . In the limit q M <C fc M the main contribution to the numerator and 
denominator comes from the first term, 

gA, 2(^/fc a/ ) 2 [s(l)a(2)-a(0)a(3)] 2 [a(l)a(2) - a(0)a(3)] 
h ~ [ S 2 (2)- S (0) S (4)] ~ [ S =(2)- S (0) S (4)] ' 

reproducing Eq. ([T4|) . 



IV. NUMERICAL SIMULATIONS 

The ?7,ew phase time formula (|17p has been tested for an incoming gaussian wave packet, 



exp[-(/c — fc ) d /4] for < fc < fc AJ , 
otherwise , 



with a localization d = 10/fc M , with a momentum distribution centered at /c — k M /2, and with an 
upper limit for the momentum distribution given by k M = KeV /he. The incoming electrons move in 
the free region before the barrier with velocity 



v a = hk M /2m = \jE M j2m = 10 c. 
For a potential barrier of height V = E M (= leV), the transit velocity is given, see Eq. ([T5|) . by 

-3 

V tra =4.5-10 C . 

Numerical data are presented in Table 1. The time in which the transmitted peak appears in the free 
region after the barrier is calculated for different values of L. For increasing L, the transit velocity, 
v tra , tends to a constant value which is in excellent agreement with the analytic value obtained from 
Eq.G5). 

To complete our numerical analysis, we have calculated the transit velocity as a function of k M L 
for different ratios of V /E M . The plots in Fig 1 clearly show that such a velocity tends to a constant 
value for V — > E M . The standard phase time, the new phase time and the numerical data are plotted 
in Fig. 2. The new phase time is in excellent agreement with the numerical simulations. The standard 
phase time represents a good approximation for increasing values of V /E M . 
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V. CONCLUSIONS 



The growing interest in understanding tunneling times in quantum physics stimulated the authors in 
looking for a new analytic formula of the phase time for wave packets transmission through opaque 
barriers. After a brief review of the derivation of the phase time formula, which is based on the 
stationary phase method, we discuss some intriguing questions about its appropriate use. In the 
opaque limit, the filter effect is responsible for a shift of the mean value of the transmitted momentum. 
This allows to compute directly the transmitted electronic density, and, consequently, by taking the 
time derivative of this density, to find the time in which the wave packet appears in the free region 
after the barrier potential. The new formula for the phase time is in excellent agreement with the 
numerical simulations and clearly shows in which cases the standard phase time, calculated by the 
stationary phase method, represents a good approximation for the transit time. The most important 
goal of the paper is the proof that, for wave packets whose upper limit of the momentum distribution 
is very close to the barrier height, the phase time is proportional to the barrier width. 

Finally, we hope that this work will find readers not only among the physicists interested in 
tunneling phenomena but also among the specialists in related branches of natural sciences which use 
the stationary phase method in their practical research. 
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L 


t 

[h/v ] 


Vtra 


Ana/Num 
o 


[h/V2mV ] 


yv /2m] 


50 


10.20 


4.9013 


91.81 % 


100 


21.41 


4.6715 


96.33 % 


150 


32.37 


4.6338 


97.11 % 


200 


43.28 


4.6209 


97.38 % 


250 


54.17 


4.6150 


97.51 % 


300 


65.05 


4.6118 


97.58 % 


350 


75.92 


4.6099 


97.62 % 


400 


86.79 


4.6086 


97.64 % 


450 


97.66 


4.6078 


97.66 % 


500 


108.53 


4.6072 


97.67 % 



h/ V2 m V « 2 A 

H/V w 0.66 • 10- 15 sec 
y/V /2m w 10^" 3 c 



Table 1: Numerical data for a localized (d = 20 A) non relativistic (v = 10~ 3 c) electron with E M = V a 
which tunnels across a barrier potential, V = 1 eV. The time in which the transmitted peak appears 
in the free region after the barrier (second column), is calculated for different values of L. It is clear 
the linear dependence on L of the transmission time. Consequently, the barrier transit velocity (third 
column) defined by L/t, tends to a constant value. The ratio between analytical and numerical transit 
velocities is given in the last column, and shows, as expected, an excellent agreement for increasing 
values of L. 
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Fig. 1 
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Figure 1: Transit velocities for an incoming (electron) wave packet centered in k = k M /2 with 
localization determined by k M d = 10 as a function of k M L. For V — > E M , the transit velocity tends 
to a constant value. 
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Figure 2: The new phase time formula (New), the standard phase time formula obtained by using the 
stationary phase method (Spm) and the numerical data (Num) are plotted as functions of y/V /E M 
for a barrier width determined by k M L = 100. The new phase time formula is in excellent agreement 
with the numerical analysis. For increasing V /E M , the standard phase time formula represents a 
good approximation. 
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